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ABSTRACT 


The  influence  of  orthotropicity  and  internal  pressure  on 
the  natural  frequencies  of  a  thin  cylindrical  shell  whose  ends 
are  freely  supported  is  found  using  the  shell  equations  of 
Washizu.  Other  common  shell  theories  are  also  presented. 

Calculations  are  performed  for  four  constant  thickness 
orthotropic  cylinders,  two  stiffened  cylinders,  and  a  basic 
isotropic  "cylinder.  For  these  calculations,  only  the  special 
cases  of  axial  and  circumferential  stiffening  were  considered. 

It  is  noted  that  there  is  a  difference  in  the  behavior  of 
the  frequencies  for  the  n  *  0,  n  *  1,  and  n  ^  2  modes.  For  the 
n  ;>  2  modes,  the  lowest  frequency  is  associated  with  predomi¬ 
nantly  radial  motion.  For  the  n  *  0  mode,  the  lowest  frequency 
is  either  associated  with  uncoupled  torsional  motion  or  pre¬ 
dominantly  radial  motion,  depending  on  the  axial  wave  length. 

For  the  n  =  1  mode,  the  lowest  frequency  is  associated  with  a 
combination  of  circumferential  and  radial  motion. 

For  the  n  2  and  n  *  0  modes,  circumferential  stiffening 
generally  increased  significantly  the  frequency  of  the  pre¬ 
dominantly  radial  mode,  while  axial  stiffening  had  little  ef¬ 
fect  on  these  frequencies.  For  the  n  *  1  beam-type  mode,  axial 
stiffening  increased  the  frequency,  while  circumferential 
stiffening  has  little  effect;  on  the  frequency. 

Internal  pressure  significantly  increased  the  frequency 
of  the  predominantly  radial  mode  for  all  cases  of  stiffening 
for  n  ^  2,  while  there  was  little  effect  on  the  n  *  0  mode. 
Internal  pressure  had  small  effect  on  the  n  *  1  beam- type  mode 
for  all  cases  of  stiffening. 
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1.  INTRODUCTION 


Many  papers  have  been  written  on  the  free  vibrations  of  a 
thin  circular  cylindrical  shell.  For  example,  Arnold  and 
Warburton  (Ref.  1)  have  investigated  the  frequency  spectrum  for 
freely-supported  circular  cylindrical  shells  using  a  set  of  dis¬ 
placement  equations  of  motion  which  are  similar  to  those  of 
Goldenveizer  (Ref.  2). 

The  effect  of  initial  circumferential  and  axial  stresses 
on  the  vibration  frequencies  of  freely-supported  isotropic  cylin¬ 
drical  shells  has  been  studied  by  Reissner  (Refs.  3,  4).  In 
Ref.  3, he  derived  a  simple  expression  for  the  frequency  on  the 
basis  of  Maguerre's  shallow-shell  theory  by  omitting  the  effect 
of  circumferential  and  longitudinal  inertia  terms.  The  results 
obtained  in  Ref.  3  cannot  be  expected  to  be  valid  for  modes 
having  a  small  circumferential  wave  number,  due  to  the  shallow- 
shell  approximation.  In  Ref.  4,  Reissner  developed  the  frequency 
equation  for  isotropic  cylindrical  shells  using  membrane  theory, 
including  all  inertia  effects. 

Fung,  Sechler  and  Kaplan  (Ref.  5)  have  also  studied  the 
effect  of  pressure  on  the  vibration  frequencies  of  freely- 
supported  cylindrical  shells.  To  include  the  effect  of  pressure, 
the  nonlinear  equations  of  Timoshenko  (Ref.  6)  are  used.  How¬ 
ever,  the  calculations  were  carried  out  using  the  frequency 
equation  of  Reissner  (Ref.  3). 

More  recently,  Armen&kas  (Ref.  7)  has  investigated  the 
effect  of  initial  stress  on  simply -supported  isotropic  cylindrical 
shells  using  a  shell  theory  developed  by  Herrmann  and  Armen&kas 
(Ref.  8)  which  includes  the  effect  on  the  motion  of  the  change 
due  to  deformation  of  the  magnitude  and  direction  of  the  applied 
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initial  io«3. 

This  investigation  is  primarily  concerned  with  the  effect 
of  circumferential  and  axial  stiffening,  as  well  as  internal 
pressure  on  the  frequencies  and  mode  shapes,  of  a  freely-supported 
cylinder. 
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2.  BASIC  SHELL  EQUATIONS  AND  RELATIONS 


Shown  in  Fig.  la  is  the  middle  surface  of  an  orthotropic 
cylindrical  shell.  The  x-y  axes  are  the  geometric  axes,  while 
the  x'-y'  are  the  principal  axes  rotated  an  angle  cp  with  respect 
to  the  geometric  axes.  The  stiffened  axis  of  the  cylinder  will 
be  taken  as  x' .  The  stress  resultants  and  couples  are  shown  in 
Fig.  lb. 

The  boundary  conditions  used  here  are  those  called  freely- 
supported  by  Arnold  and  Warburton  (Ref.  1)  and  are  given  by 

=  V(J>  0'*)  =  ° 

W-(0j4,i)  i)  —  O 

(1) 

M =  M„/7,  t>,t)  =  o 


The  effect  of  other  boundary  conditions  on  the  modal  character¬ 
istics  of  thin  isotropic  cylindrical  shells  has  been  investigated 
in  detail  by  Forsberg  (Ref.  9). 

2.1  Thin  Cylindrical  Shell  Equations  using  Washizu's 
Shell  Theory 

From  Fig.  1,  assuming  that  the  thickness  is  small  compared 
to  the  radius  (we  limit  ourselves  here  to  Love's  first  approxima¬ 
tion,  i.e.,  1  +  z/R  ~  1),  the  stress  resultants  are 
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(2) 


(3) 


and  the  stress  couples  by 


(4) 


(5) 


(6) 


(7) 
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Strain-Displacement  Relations 

The  effect  of  transverse  shear  is  neglected  and  normal 
strains  are  assumed  to  be  zero.  When  the  Euler-Bemoulli-Navier 
hypothesis  is  also  employed,  the  strains  can  be  written  as 

4c  *  Sr.-fcfr*  (8) 

€&-  ~  ~  *  ^9  (9) 


~  Ken  -  A*  h] 


xfr 


(10) 


where  t  and  are  the  mid -surface  strains  and 

Kx  ;  tie-,  and  are  the  curvature  strains. 

Using  what  Washizu  (Ref.  10)  calls  the  thin  shell  of 
Type  D,  the  strains  are  related  to  the  mid -surface  displacements 
by 


~h 


± 

Z 


(11) 


z  /  ?iT  ur  ^  / 

^  -  *  » '  1 


f2l 


(12) 


YX)  -  21  +  L  in  +  ± 

'*  9x  A  9*  R 


>>  +  Ilf  21-  ur)+2l[2i+i 

9x  9$  0x(  J  yk  J 


(13) 
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(14) 


9V 


( 


Tur  i  2JL\ 


(15) 


J_  ( £jr  2* 
K  \dx2fi  gx 


'■) 


(16) 


These  expressions  include  nonlinear  terms  in  midp  lane  strains, 
but  not  in  the  curvature  strains  e* 

Stress-Strain  Relations 

The  stress -strain  relations  for  a  thin  orthotropic  cylin¬ 
drical  shell  are 


Nx  ~  £k9  +C:z  €&0  +  3  fa+o 

(17) 

N*  —  ^/z.  £xa  -+  Czi  6#o  3-  CZ5yx»0 

(18) 

N**  -  Cl3  +  ^i-3  +  ^33  fa&c 

(19) 

Mx  -  -  (  Pn  Kx  +  T>iz  K»  +ZP,3  fr*#) 

(20) 

Mi  -  ~~  { P/2.  fSx  +  Pm-  Kxi) 

(21) 

-* 

The  nonlinear  terms  in  the  curvature  strains  are  omitted  here 
since  they  will  offer  no  contribution  to  the  final  linearised 
cylindrical  shell  equations  under  initial  membrane  stresses. 
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Mfp  —  O'V?  +  -^23  ^-^33 


(22) 


Equations  of  Equilibrium 

For  the  special  case  of  the  boundary  conditions  given  in 
Eqs.  (1),  the  principle  of  virtual  work  as  applied  to  the  cylin¬ 
drical  shell  can  he  written  as 


Nxo6xo  ~t  N&-  S  6»c 


+  NxtSh i* 


Mx  <f  fa  ~ 


~  -Z.  Mx&  &  fax*  jRjxJ#- 


-h  f^Sif 


-tf=±  ur)  R  JxJ&  -  o  (23) 

where  S  is  the  area  of  the  raid-surface  and  Fx,  F0>  and  Fz  are 
the  body  forces  per  unit  of  mid -surface  area.  The  equilibrium 
equations  obtained  from  Eq.  (23)  depend  on  the  strain-displace¬ 
ment  relations  given  by  Eqs.  (11)  through  (16).  This  leads  to 
the  following  nonlinear  equilibrium  equations  of  Washizu  (Refs. 
10  and  11) . 


4  FXR  =  o 


(24) 


7 


I[N.RI ¥  + 

if-  i)+ 

+ N»f]  -  f”’)  <2» 


*  *4»^  89] 


*^]+  “'(!**&-$)  +  ""ft 


-j__2.  [>  QMx  +~  ®My*-  -x  J?  _i. (- 

a*  l  0*  <9*  J  m/e  / 


*  j?  r_i.  Qm*-  .f 


mx?  o+ 


(26) 


The.  Frequence 


mat  ion 


For  the  vibration  problem,  the  nonlinear  equilibrium  equa¬ 
tions  [Eqs.  (24)  through  (26)]  are  used  as  a  basis  for  deriving 
the  influence  of  internal  pressure  on  the  frequencies  and  mode 
shapes  of  the  thin  cylindrical  shell.  The  shell  stresses  and 
displacements  can  be  expressed  in  terms  of  deviations  from  an 
initially  stressed  configuration.  We  take  here  the  case  where 
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the  only  initial  stress  resultants  are  ^  and  and  all  the 
initial  stress  couples  are  zero.  That  is,  let 

— ■— 

N*  -  N&  +  (27) 


Nxa  =  /Vjra 


A'/x  "  Afx  ,  / 


(28) 


and 


^  =  u  +  u 


'V~-  V  +  4T 


ur w  ur  4-  ur 


(29) 


where  Nx,  N0,  u,  v,  and  w  are  those  stress  resultants  and  dis¬ 
placements  resulting  from  the  static  loading  conditions  and 
8X»  Nx0,  H0,  u,  v,  and  w  are  the  vibratory  stress 

resultants,  couples,  and  displacements. 


Considered  in  this  report  is  the  case  in  which  the  internal 
pressure  is  taken  as  Constant  Directional,  i.e.,  the  direction  and 
magnitude  per  unit  original  area  remains  unchanged  during  deforma¬ 
tion.  For  a  consideration  of  the  case  where  the  direction  and 
magnitude  per  unit  original  area  changes  with  deformation  due  to 
the  internal  pressure,  see  Ref.  8.* 


All  the  vibratory  stress  resultants  and  couples,  as  well 
as  the  vibratory  displacements,  are  considered  infinitesimal 
quantities  for  the  vibration  problem  and  it  is  justifiable  to 
neglect  their  squares  and  products.  It  is  also  assumed  that  Hx 
and  Ng  are  constant  and 


*For  that  case,  additional  terms  p  ^  (^  +  v),  and 
~  p[|~  +  ^  (|^  —  j|)]  would  be  added  to  the  left-hand  sides 
of  Eqs .  (31),  (32),  and  (33),  respectively. 
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/V1 


£V 


F  = 


,2V 

*■ 


(30) 


With  this  simplification  Eqs.  (24)  through  (26)  become  (here,  for 
convenience  of  notation  we  now  write  the  vibratory  stresses  and 

O  A/  A* 

displacements  Nx,  N0,  u,  v,  and  w,  as  Nx>  N0,  u,  v,  and  w). 


'2%  t  L  1***  +  ur  +  =  o 

9x  z  9*  0r‘  Kl  0*' 


(31) 


N± 

R 


V'Mx  +  * 
2xx  /? 


,  A /♦  /  2<r  . 

+T‘(pf 


/ 

+  — 

0M& 

+  N*  3‘f 

0&1- 

0x*- 

'  9jr_ 

-  ur) 

>  0  1  u/~ 
—  W* 

96- 

J 

=•  <9 


(33) 


where  Eqs.  (31)  through  (33)  represent  the  equilibrium  equations 
of  free  vibration  of  a  thin  cylindrical  shell  about  an  initially 
stressed  configuration. 

By  combining  Eqs.  (11)  through  (16)  and  Eqs.  (17)  through 
(22)  and  substituting  this  result  into  Eqs.  (31)  through  (33),  the 
following  differential  equations  in  u,  v,  and  w  are  obtained. 
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C„  T<<  a  Qs  -3JJL  ,  ±  £n  Tu  +  ?1£ 

Z'.K'  «  cix  0x9»  **  ^  ^  *** 

/  £/*  -#■  £j/  p  V"  _J_  _/  ^ z  j  £  **r _ f_  Crz  QjC 

+  y?  <T«  ■  9x9*  C,C  0ex  X  9x 

__/  Ctt  tj£  +  _;_  Pxx  FJU.  V)*r)  4-  G,(fi»,  N*,",  r,  “) 

/?*  ai  if-  fi  cu 

M*  QlC1  _  (34) 

cU  9fi 


\  c„  7iu  .  /  Pi  +  csi  ?y  ,  ±  £«  iLk:  _,  £" 

fix  9xl+  X  ai  xx  cA  0ex  cU  #xl 

z  1^3  2  V  .  /  £V  _  1  Qi.  £J£-  ±  £*f 

4  x  fix  9x9»  fi  Cxi  ?»"  X‘  Cxi  2x  *'  ci-  9* 

+  _L  =2”  Fx  (X,  v,  nr)  +  Gx( fix,  fi+,  ",  ^nr) 

K‘  cU 

_  *fi  2-f  =  n  (35) 

fix  M 


ii 


J_  &Z  Py  J if  /  )  V  <^2i.  '&  if 

4  Cll  9x  '  /f "  ^  yf  <^'  'A  /f*  HI  A 


/  <^2Z  _  _/  .A  2. 

X'  <r*i  ,f‘  c2: 


P/t  j?'-  PV  .  n  fi/x.  +-  A  Pi  i  p^4 /■ 

J>«  A«  j>'  £*»*' 


+4, 


'zz 


Pn/~ 

Pu  p^ 


A 

<f,I 


f^r.  ^  Fi  (">''>“)  +  6,  (W*,M> 


(36) 


t>£  9V 


£> 


The  expressions  for  Fp  F2,  F^,  and  Gp  G2,  G^  appearing 
in  Eqs.  (34)  through  (36)  are  given  in  Table  1.  The  explicit 
values  of  Fp  F2,  F^  used  here  are  those  of  Washizu  (Ref.  10), 
Novozhilov  (Ref.  16),  and  Goldenviezer  (Ref.  2),  while  the  Gp 
G2,  are  those  of  Washizu  (Ref.  10)  and  Herrmann  (Ref.  8). 

Numerous  other  versions  of  the  cylindrical  shell  equations 
have  been  published  in  the  literature.  All  of  these  versions  may 
be  written  in  the  form  of  Eqs.  (34)  through  (36)  in  which  minor 
differences  are  reflected  in  the  functions  Fp  F2,  F^,  and  Gp 
^2*  ^3*  Some  other  common  theories  are  also  listed  in  Table  1. 

We  will  be  concerned  in  what  follows  with  the  special 
cases  of  orthotropicity  where  the  geometric  and  principal  axes 
are  parallel  and  perpendicular,  i.e.,  cp  »  0  and  90  degrees.  This 
corresponds  to  the  conditions 


<^/3 

cU 


63  _  A?  _  2>is  __  q 

Cxr.  A't  -Al 


(37) 


With  this  simplification,  a  solution  to  Eqs.  (34)  through 
(36)  which  satisfies  the  boundary  conditions  [Eq.  (1)]  is 
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(38) 


n  =  Um„  Cos  *ZJp<  l 
flT  =  /#»*  si*£Z£ 

W-  urm,  si»2Ll£. 

where  m  and  n  are  integers  and  w  is  the  natural  frequency  of 
vibration.  The  integer  n  represents  the  number  of  circumferential 
waves,  while  the  integer  m,  the  number  of  axial  half  waves.  A 
pair  of  integers  m,  n  then  specifies  a  particular  mode  shape. 

Upon  substituting  Eqs.  (38)  into  Eqs.  (34)  through  (36) 

there  are  obtained  three  homogeneous  linear  algebraic  equations 

in  u _ ,  v.  and  w  .  the  nontrivial  solution  of  which  requires 

mn  mn  mn 
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-/•  /TV  (tf  *•+  0  “  %/l  *  ( 39 ) 


where 


f\  -  tfojwvfcil 


JVI* 


n'JTA! 

J. 


(40) 


Note  that  the  vibration  frequency  and  initial  stresses  are  ref¬ 
erenced  with  respect  to  C^*  i-e-f  the  stretching  stiffness 
along  the  unstiffened  y'  direction. 

When  n  and  vanish  and  the  material  is  isotropic, 

Eq.  (39)  reduces  to  that  given  by  Arnold  and  Warburton  (Ref.  1). 
If  the  determinant  is  written  as 


-  yr 

#3  I 


OLn, 

Ou-L-ft 


a 


13 


a 


3Z 


#23 

#33~  fC 


-  o 


there  is  provided  a  frequency  equation  of  the  form 


(41) 


(42 


where 


K  —  +#22  +  #33 


K 


0  1! 

(X  n 

Cl  13 

cizi. 

#23 

— 

4- 

-V- 

#2-f  (XVL 

Ctj< 

#53 

#32 

#33 

(43) 
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dll 

(Ln. 

^13 

a2, 

CLix 

^.3 

d-3! 

(X^z. 

&3J 

Associated  with  the  three  positive  roots  (ft^  <  <  ft-) 

of  Eq.  (42)  are  three  inode  shapes  given  by 


where 


^">1 


irimij 


J>  = 

N,  . 

Nv  - 
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(44) 


(45) 


The  general  character  of  these  mode  shapes  depends  on  the 
shell  parameters  and  nodal  pattern. 
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3.  APPLICATION  TO  CONSTANT  THICKNESS  ORTHOTROPIC  AND 


STIFF 


CYLINDERS 


The  results  of  the  preceding  section  are  now  applied  to 
an  orthotropic  cylinder  of  constant  thickness  and  also  to  a 
stiffened  cylinder  of  skin  thickness  t  ,  which  is  integrally 
stiffened  on  one  side  (Table  4).  In  the  latter  case,  it  is 
assumed  that  the  spacing  of  the  stiffeners  is  small  compared  to 
the  radius  and  length,  so  that  for  the  purpose  of  studying  over¬ 
all  or  average  behavior,  it  is  plausible  to  replace  the  actual 
cylinder  by  an  equivalent  uniform  orthotropic  cylinder.  Note 
that  in  this  report  the  shell  theory  used  applies  to  cross  sections 
of  shell  elements  which  are  symmetrical  with  respect  to  the  middle 
surface.  For  asymmetric  cross  sections,  there  will  be  coupling 
between  the  mid-surface  strains  and  changes  of  curvatures.  This 
coupling  effect  is  omitted  here. 


3.1  Properties  of  Constant  Thickness  Orthotropic 
Cylindrical  Shell 

For  an  orthotropic  material  of  constant  thickness  h,  it  is 
possible  to  express  the  stress-strain  law  in  the  form 

( T\  —  b„  6x  +  &/*■  £&  + 

OV  -  b/i.  6x  +  bzi  6&  -f  (46) 

&X&  =  4,3  +  4,  ~h  &  x  3 


where  the  b. .  are  the  elastic  stiffnesses  of  the  cylindrical 
shell  referred  to  the  x-y  axes.  The  stiffness  b^j  can  be  written 
in  terms  of  the  stiffnesses  referred  to  the  principal  axes, 
x'-y’,  as  shown  in  Ref.  12  as  follows: 
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bff  —  b„  Cos*#  +  bzi  St/t*/  -f-  (&SZ  +  A &33J  Cos*/ s"?*  / 
b/i  ~  (bn  +  ^2V  ~  1  b33^  Cos*  /  S/J9  V  +  £/z  fc*S  V  -+S"7*<P J 


(bf(  +&21  —  3  4/2.J  Cos  */  s/'t1/’  £33  (Cos  L/~  ">/*/*/> J 


where  the  £our  independent  material  constants  b!.  have  the  foil ow- 
ing  physical  meaning 


The  stiffnesses  of  the  cylinder  are  completely  determined  by 
the  moduli  of  elasticity  E^,  the  shear  modulus  Gx»g»,  the 
Poisson's  ratio  v0»x«,  and  by  the  angle  9. 

When  the  geometric  and  principal  axes  are  parallel,  i.e., 
9=0  deg, 
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(49) 


When  the  geometric  and  principal  axes  are  perpendicular, 
i.e. ,  9  *  90  deg 

-  izi  j  ^/z  =  /  &zz  ~ 

(50) 

^33  ~  &33  j  &/3  ~  ^2  3  ~  & 

When  the  material  is  isotropic. 


?)£ 


£ 

M'+v) 


(51) 


Note  that  for  the  constant  thickness  orthotropic  case,  the 
stretching  stiffnesses  and  the  bending  stiffnesses  are  not  in¬ 
dependent,  i.e.,  from  Eqs.  (2)  through  (7)  and  (17)  through  (22), 
they  are  given  by, 

Dij  -  ~  k/  h  Cij  —  k/  A  (52) 

3.2  Properties  of  Stiffened  Cylinder 

A  cross  section  of  a  repeated  element  of  the  stiffened 
cylindrical  shell  is  shown  in  Table  4.  The  stiffeners  are 
aligned  along  the  x’  axis  and  make  an  angle  9  with  respect  to 
the  x  axis. 


The  stretching  stiffnesses,  referred  to  the  principal 
axes  x'-y',  are  C^,  and  may  be  approximated  (with 

the  assumption  of  no  stress  lag)  by 


C/Z  =  £1/  =  7?  Cj 


[±t_t 

I  -f  f  f  ~  ] 

6X  6i\/  +  £3/6* J 


(53) 


<TZZ  =  Cs 


I  -L_  btv 

‘  Ay  is 


j  -j-  —  —  f 

^s\!+  /iw 


Sy 


<?*s 


it 

I  \  ^  f iw/&y\ 

h  /+  ti  //*/ J 


where  Ce  is  the  stretching  stiffness  of  the  basic  unstiffened 
s 

skin 


4  =  (54) 

The  bending  stiffnesses,  referred  to  the  principal  axes 
x'-y',  are  D^,  Dj^*  D^,  D33*  an<*  b®  aPProx^-m£te^  (see 
Refs .  13  and  14)  by 
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D,*  =  K 


4  = 


_ 


/z 


a_  C  _  .  J>s  (/-v) 

JJ>y 


A 


\-h  £  —  f±* 
1  *  ( *, 


(55) 


where  C  is  the  torsional  rigidity  of  the  web,  P  is  a  constant 

based  on  b  and  t.  which  varies  from  0.333  to  0.141,  and  De  is 
w  w  s 

the  bending  stiffness  of  the  unstiffened  skin 


/~2?z 


(56) 


The  thin  cylindrical  shell  with  closely-spaced  stiffeners  will  be 
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considered  as  an  equivalent  orthotropic  thin  cylinder  with  given 
independent  bending  and  stretching  stiffnesses.  With  this  ap¬ 
proach  and  in  conjunction  with  the  stiffness  transformations 
given  by  Eqs.  (47),  the  stiffnesses  C.j  and  referred  to  the 
geometrical  axes  x-y  can  be  written  in  terms  of  the  stiffnesses 
C!.  and  D!.  referred  to  the  principal  axes  x^y*  as 


(Cn j  X>„\  ~  (^n  ,  D/t'jZ-os  V  +  (&z  2  /  J  5/jt 

■+■  2.  (  tj  \ +  X  (  ^33y  $33 ()  CoS  V  jS/''7  *0 

(fl*/  ^/2^j  ~  (  Cn,  J>//J  +  (<?2l,  Pzi)  ~  *  (£33 j  ^33 )  C05*  & 
~t  (  C,z,  D,'i  )(C05  y  +  5/* 

(^22 j  &Z “  (  £// /  Ptf  )  3//7  *</>  +■  ( Czzj  COS*0 

+*[(  c,z,  P,t)  +Z(C&,J>*j[  Cosy  s,*v 
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(c,3f  d/3 )  =  [( v  -  f  j>»)  °°*  v 


(57) 


—  ^'2-/  Jlz  J  ~t~  *2  (£33  j  ^*3  Jj-  Cos*fl>  -  s/n  *^)\CoZ><fj  S/*t$ 


(c.1),  %)=  ( c«',  &  i'yW  <*  -  ^ ;  r:)  s,* v 


_  j  (Cn,j>,0  +  a( c,l,  ft'jj(cflsv-^V/) 


(^3Jy  —  C  Cr/  /  Ptf  )  1~  (  Cz  l  y  -A^ )  ^  ^  1 '  ) 


S/*7  V 


+  (  ("53/  T>s?)  (Q.O$X'<f>  -  S/*  yj) 


When  the  geometric  and  principal  axes  are  parallel,  i.e., 
cp  =  0  deg 
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.e. 


> 


When  the  geometric  and  principal  axes  are  perpendicular. 


“  90  deg 
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3.3  Numerical  Application 

a)  Constant  Thickness  Orthotropic  Cylinder 

Four  examp le  cylinders  were  selected  for  numerical  applica 
tion  for  the  constant  thickness  orthotropic  case  and  their  prop¬ 
erties  are  given  in  Table  1.  Calculations  for  the  isotropic  case 
were  also  carried  out  for  purposes  of  comparison.  The  thickness - 
to-radius  parameter  was  chosen  as  0.001  for  these  calculations. 
Note  that  for  the  constant  thickness  orthotropic  case,  the  ratio 
of  stretching  stiffnesses  and  bending  stiffnesses  are  the  same, 
i.e. , 


Cij  _  J)t/  _ 

Cz'l  bzi 


(60) 


b)  Stiffened  Cylinder 

Two  examples  of  stiffened  cylinders  were  selected  whose 
bending  stiffness  ratio  is  the  same  as  those  of  cylinders 

(1)  and  (5)  for  the  constant  thickness  orthotropic  case.  In 
order  to  compare  the  frequency  spectrum  for  the  constant  thickness 
orthotropic  cylinder  and  the  stiffened  cylinder  in  a  meaningful 
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way,  the  term  (l/R^XD^/CJ^)  for  the  stiffened  cylinder  is  set 
equal  to  (l/12)(h/fc)2. 
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►  .  DISCUSSION  OF  CALCULATIONS 


4.1  Constant  Thickness  Orthotropic  Cylinder 

The  general  character  of  the  frequency  spectrum  is  indi¬ 
cated  in  Figs.  4  through  13.  These  curves  indicate  that  the 
modes  for  which  n  <  2  have  a  different  behavior  than  those  for 
n  ^  2.  The  modal  characteristics  for  n  *  0  and  n  -  1  for  the 
isotropic  case  are  essentially  independent  of  the  bending  stiff¬ 
ness  of  the  shell  (i.e. ,  independent  of  h/R) ,  since  for  these 
cases  the  deformation  consists  mainly  of  stretching,  whereas 
for  high  n  the  frequency  characteristics  depend  more  strongly 
on  bending  stiffness.* 

Let  us  consider  in  detail  the  frequency  spectrum  for 
n  -  0.  For  this  case  one  of  the  three  roots  of  Eq.  (42)  repre¬ 
sents  the  frequency  of  the  uncoupled  torsional  mode,  while  the 
other  two  frequencies  involve  predominantly  longitudinal  and 
radial  motion.  From  Eq.  (39)  the  uncoupled  torsional  frequency 
is 


(61) 


while  the  torsional  frequency  of  a  thin -walled  circular  b&am 
according  to  St.  Venant  torsion  theory  is 


(62) 


For  short  axial  e  lengths  (large  raR/4),  the  effect  of  bending 
stiffness  (h/R)  tn  iv  become  important  for  n  ■  0  and  1.  (See 
Figs.  2  and  3) 
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The  other  two  frequencies  also  shown  in  Figs.  4  through  8  have 
as  asymp totes  the  frequency  of  axial  vibrations  of  a  bar. 


and  the  frequency  of  radial  vibrations  of  a  ring  in  plane  strain 
for  long  axial  wave  lengths. 


i  i’ll 


and  a  ring  in  plane  stress  for  short  axial  wave  lengths, 


(65) 


For  the  n  *=  0  case,  the  lowest  frequency  is  not  necessarily 
associated  with  predominantly  radial  motion.  Figures  4  through 
8  show  that  the  lowest  natural  frequency  can  be  associated  with 
either  the  predominantly  radial  mode  or  torsional  mode  depending 
on  mR/i.  The  value  of  mR//  for  which  the  lowest  natural  fre¬ 
quency  changes  from  the  predominantly  radial  mode  to  the  torsion¬ 
al  mode  depends  on  the  stiffness  ratio  bj^/b^  and  the  angle  q>. 
From  Figs.  4  through  8,  it  is  observed  that  the  torsional  fre¬ 
quency  is  slightly  affected  by  the  stiffness  ratios  4x^22 
and  the  angle  q>,  while  the  axial  frequency,  since  it  depends 
mainly  on  the  stiffness  in  the  axial  direction,  is  affected  by 
the  stiffness  ratio  b^/b^  only  when  9  »  0  deg.  It  also  may 
be  noted  that  the  rise  in  the  predominantly  radial  frequency  at 
very  short  axial  wave  lengths  (large  raR//)  is  due  to  the  effect 
of  h/R. 
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Let  us  turn  now  to  the  case  for  n  c  1.  Here,  for  long 
axial  wave  lengths  (small  mR//) ,  the  mode  associated  with  the 
lowest  frequency  involves  strong  coupling  in  the  circumferential 
and  radial  directions.  For  long  axial  wave  lengths,  this  cor¬ 
responds  to  the  vibration  of  a  cylinder  as  a  beam.  The  asymptotes 
for  cases  (1)  and  (5)  that  are  shown  in  Fig. 14  are  those  of  a 
simply-supported  thin-walled  circular  cross -sectioned  beam  as 
found  from  elementary  beam  theory  to  be 


(66) 


Figures  9  through  13  show  that  for  long  axial  wave  lengths 
(small  mR/Z)  circumferential  stiffening  has  a  negligible  effect 
on  the  beam-type  frequencies  while  axial  stiffening  has  a  signifi 
cant  effect;  whereas  for  short  axial  wave  lengths  (large  mR/l) 
axial  stiffening  has  a  slight  effect  and  circumferential  stiffen¬ 
ing  has  a  significant  effect  on  the  frequency. 

Now  let  us  consider  the  case  for  n  ^  2.  The  asymptotic 
values  of  the  three  frequencies  for  long  axial  wave  lengths 
(small  mR/i)  are  (for  nx  *=  n0  *  0): 

Flexural  vibrations  of  a  ring  (e^  *  0,  u  =  0) 


1 
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Axial  shear  vibration  (w  *  v 


0) 


A 


(68) 


Extensional  vibrations  of  a  ring  (<e  *  0,  u  **  0) 


For  the  n  2  cases,  the  lowest  frequency  corresponds  to  the  pre¬ 
dominantly  radial  mode  for  all  values  of  mR/i  and  stiffnesses, 
while  the  remaining  two  correspond  to  predominantly  axial  and 
circumferential  motion.  In  general,  these  two  are  at  least  an 
order  of  magnitude  higher  than  the  predominantly  radial  fre¬ 
quency.  From  Figs.  9  through  13,  the  effect  of  the  stiffness 
ratio  and  the  rotation  of  the  principal  axes  of  stiffness  through 
an  angle  of  90  degrees  can  be  studied.  It  may  be  seen  that 
stiffness  ratio  ^as  effect  on  the  predominantly 

radial  frequency  for  long  and  intermediate  axial  wave  lengths 
when  <p  »  0  deg  (axial  stiffening),  while  this  frequency  shows  a 
marked  increase  when  the  stiffness  ratio  bj^/b’ ^2  increases  for 
<p  «  90  deg  (circumferential  stiffening).  Specifically,  it  may 
be  noted  that  for  raR/i  “  0.01,  n  «  3,  bj^/b^  "  24.2,  when  cp  « 

0  deg  the  frequency  is  0.00225,  and  when  <p  *  90  deg,  the  fre¬ 
quency  is  0.0108;  while  for  the  corresponding  isotropic  case, 
the  frequency  is  0.00219.  For  all  values  of  stiffening  both  in 
circumferential,  as  well  as  axial  direction  for  large  axial  wave 
lengths,  the  frequency  is  proportional  to  n  (for  sufficiently 
high  n)  as  indicated  by  Eq.  (67). 
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In  all  the  cases  of  stiffening  for  n  ^  2,  the  value  of  n 
for  which  the  minimum  frequency  occurs,  increases  with  increasing 
mR/i.  This  is  also  true  for  the  isotropic  case. 

4.2  Stiffened  Cylinder 

In  the  preceding  section,  the  results  were  obtained  from 
the  analysis  of  a  cylindrical  shell  which  was  orthotropic  due  to 
the  natural  properties  of  the  material.  Here  we  will  consider  a 
shell  which  is  orthotropic  due  to  attachment  of  stiffeners  on  to 
the  basic  thin-walled  cylinder.  The  purpose  of  this  analysis  will 
be  to  compare  the  frequency  characteristics  obtained  from  a  study 
of  the  stiffened  cylinder  with  those  obtained  from  the  constant 
thickness  orthotropic  cylinder. 

Figures  15  through  18  show  the  results  of  the  stiffened 
cylinder  analysis.  The  properties  of  the  cylinders  are  given  in 
Table  4.  The  essential  difference  in  the  properties  of  the 
stiffened  cylinder  and  the  constant  thickness  orthotropic  cylinder 
is  that  in  the  former,  the  stretching  stiffness  and  bending  stiff¬ 
ness  ratio  are  quite  different  while  in  the  latter  case  they  are 
identical.  In  particular,  for  the  stiffened  cylinders  considered 
in  this  section,  when  the  bending  stiffness  ratio  D{]/D22  was 
taken  to  be  the  same  as  that  of  the  bending  stiffness  ratio  in 
Cases  1  and  5  of  the  constant  thickness  orthotropic  cylinder, 
namely  D{]/D22  *  24.2,  the  stretching  stiffness  ratio  Cll/C22  in 
the  case  of  the  stiffened  shell  was  1.26.  Returning  to  Figs.  15 
through  18  one  finds  that  for  the  n  ■  0  and  n  *  1  cases,  which 
correspond  to  predominantly  stretching  of  the  shell  middle  surface, 
the  values  of  the  frequency  are  approximately  the  same  as  those  of 
the  isotropic  case.* 

# 

For  <p  *  0  deg  and  n  =  0,  there  is  a  slight  increase  in  the  fre¬ 
quency  (when  compared  with  the  isotropic  case)  of  the  mode  associ¬ 
ated  with  predominantly  axial  motion  since  this  frequency  is  ef¬ 
fectively  proportional  to ]Cii/C22 •  This  increase  is  also  true  for 
the  n  »  1  case  when  <p  *=  0  and  mR/£  >  5,  due  to  the  slight  effect 
of  bending  for  these  high  values  of  mR /&. 
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For  the  n  ^  2  cases,  the  effect  of  bending  becomes  ex¬ 
tremely  important  even  for  n  as  low  as  2.  If  we  compare  Fig.  17 
j  with  Fig.  9  (which  correspond  to  <p  **  90  deg),  it  is  evident  that 
!  the  values  of  the  t> redominant ly  radial  frequency  for  the  stiffened 

i 

|  cylinder  are  approximately  the  same  as  those  of  the  constant 
j  thickness  orthotropic  cylinder  when  mR /i  <  0.5,  and  n  ^  2.  For 
!  higher  values  of  mR/i,  the  frequency  of  the  stiffened  cylinder 
I  decreases  below  that  of  the  corresponding  constant  thickness 

I 

orthotropic  cylinder  for  all  values  of  n  ^  2.  Note  that  this 
]  decrease  diminishes  with  increasing  n,  so  that  for  very  high  n, 

|  the  frequencies  for  both  thesr  cylinders  (constant  thickness  and 
stiffened)  again  become  approximately  the  same.  This  is  ex- 
I  plained  by  the  fact  that  for  large  values  of  n  and  mR/£  the  in- 
|  fluence  of  bending  is  predominant.  When  the  stiffeners  are 
I  aligned  in  the  axial  direction  (cp  -  0  deg),  one  observes  from 
s  Figs.  13  and  18  that  for  n  ;>  4,  frequencies  for  both  type  cylinders 
•  are  nearly  the  same  for  long  axial  wave  lengths  (low  mR/i);  for  in¬ 
termediate  axial  wave  lengths,  the  difference  in  the  frequency  for 
these  cylinders  becomes  appreciable,  while  for  short  axial  wave 
lengths,  this  difference  again  becomes  small.  For  n  =  2  and  3, 
the  frequency  of  the  cylinder  with  axial  stiffeners  is  lower 
than  the  corresponding  constant  thickness  orthotropic  cylinder 
for  all  but  short  axial  wave  lengths. 

4.3  Effect  of  Internal  Pressure 

For  a  closed-end  cylinder  under  internal  pressure  p, 

n*  =  /  fix  =  (70) 

In  order  to  gain  some  physical  insight  into  the  magnitude  of  the 
pressure  stiffness  parameter  n0,  one  can  express  ne  in  terms  of 
the  initial  static  circumferential  stress  as  follows, 
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/ 


(71) 


(T, 


For  the  constant  thickness  orthotropic  cylinder 


E* 


y>* 


With  vgtx i  *  vx'e»  =  0*3,  Eg i  =  10^  psi.  This  becomes 

0$  ~  A/  XIO1  tfe  fsc\ 

Shown  in  Figs.  19  through  21  is  the  variation  of  with 
the  pressure  stiffness  parameters  n0,  nx  for  the  isotropic  case* 
for  three  values  of  the  axial  wave  length.  Comparison  of  these 
figures  shows  that  the  value  of  n  for  which  the  lowest  frequency 
begins  to  vary  significantly  with  n0  depends  on  the  axial  wave 
length  wR/H.  For  mR/i  =0.06  there  is  a  significant;  increase  in 
with  n q  when  n  ^  2,  for  mR/^  =  0.5  the  increase  becomes  signifi¬ 
cant  when  n  ;>  5fand  for  mR/i!  =  3  when  n  ^  10.  For  all  these 
values  of  mR/i,  the  two  larger  frequencies,  which  correspond  to 
predominantly  tangential  motions,  are  little  affected  by  n0,  nx» 

a)  Constant  Thickness  Orthotropic  Cylinder 

The  effect  ci:  internal  pressure  on  the  frequencies  of  the 
axisymmetric  mode  (n  =  0)  is  illustrated  in  Figs.  4  through  8. 

For  both  <p  =  0  and  90  degrees  and  all  values  of  stiffening,  the 
frequency  of  the  predominantly  radial  mode  is  slightly  increased 
by  the  addition  of  internal  pressure  at  high  values  of  mR/iA  The 
increase  in  the  frequency  of  the  torsional  mode,  due  to  pressure, 

*In  all  the  pressure  investigations,  nx  is  always  included,  and 
is  taken  equal  to  tig/2. 


32 


increases  with  increasing  b^/b^  and  is  the  same  for  cp  «  0  and 
90  degrees.  This  increase  in  frequency  diminishes  rapidly  with 
increasing  values  of  mR/i.  Internal  pressure  has  a  negligible 
effect  on  the  predominantly  axial  frequency  when  cp  «  0  and  90 
degrees  for  all  values  of  the  stiffness  ratio 

For  the  n  **  1  mode,  the  effect  of  pressure  on  the  lowest 
frequency  is  shown  in  Figs.  22  through  26.  For  stiffening  in  the 
circumferential  direction  (q>  *  90  deg)  the  frequency  is  independent 
of  the  stiffness  ratio  ^{x/^22*  There  is  an  increase  in  frequency 
due  to  pressure  in  these  cases,  which  is  large  for  small  values 
of  mR/i  and  diminishes  with  increasing  mR/i.  For  stiffening  in 
the  axial  direction  (cp  =  0  deg),  the  frequency  does  depend  on  the 
stiffness  ratio  bj^/bJ^*  There  is  an  increase  in  frequency  in 
these  cases,  which  is  small  for  low  values  of  mR/-0,  diminishes 
for  intermediate  values  of  mR/iJ,  and  increases  for  high  values 
of  mR/A.  The  effect  of  pressure  on  the  two  higher  frequencies 
for  n  a  1  is  negligible. 

For  n  ;>  2,  Figs.  22  through  26  show  that  the  lowest  fre¬ 
quency  increases  significantly  with  internal  pressure  for  all 
cases  of  stiffening.  The  increase  in  the  frequency  generally 
diminishes  with  increasing  mR/i.  It  is  observed  that  the  in¬ 
crease  in  frequency  due  to  pressure  is  greater  for  stiffening  in 
the  axial  direction  than  in  the  circumferential  direction.  Thus 
for  mR /l  *  0,01,  bj^/b^  “  24.2,  and  n  *  2,  for  stiffening  in 
the  circumferential  direction  (<p  *  90  deg),  the  unpressurized 
frequency  **  0*00383,  while  in  the  corresponding  case  of  axial 
stiffening  (cp  ■  0  deg),  ■  0.00130.  However,  with  the  addition 
of  internal  pressure  (n0  »  2nx  *  0.001),  when  qp  ■  90  deg  the  fre¬ 
quency  **  0.0426,  while  for  cp  *»  0  deg,  “  0.0424.  This  same 
pattern  of  behavior,  namely  the  influence  of  internal  pressure 
on  the  frequency  being  greater  for  axial  stiffening  than  for  cir¬ 
cumferential  stiffening,  is  found  for  all  n  >  2  as  well  as  n  *  2. 
Also  illustrated  in  Figs.  22  through  26  is  the  fact  that  for 
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9  38  90  deg,  the  increase  in  frequency  due  to  internal  pressure 
decreases  with  increasing  values  of  the  stiffness  ratio  bil/»22> 
while  for  q>  =  0  deg,  this  increase  is  independent  of  b^j/b^* 

For  n  ^  2,  internal  pressure  has  a  negligible  effect  on  the  two 
larger  frequencies. 

b)  Stiffened  Cylinder 

The  effect  of  pressure  on  the  n  =  0  mode  for  the  stiffened 
cylinder  is  similar  to  that  of  the  isotropic  case.  This  is  due 
to  the  fact  that  the  stretching  stiffness  ratio  Cji/C^  of  the 
stiffened  cylinder  and  the  isotropic  cylinder  are  nearly  the  same. 
For  both  9=0  and  90  degrees,  the  frequency  of  the  predominantly 
radial  mode  is  slightly  increased  by  internal  pressure  for  high 
wR/i.  The  increase  in  the  frequency  of  the  torsional  mode  due  to 
internal  pressure  is  the  same  for  9=0  and  90  degrees  and  again 
diminishes  rapidly  with  increasing  mR/iL  The  axial  frequency  is 
unaffected  by  internal  pressure  for  both  9=0  and  90  degrees. 

Turning  to  the  n  =  1  case,  Figs.  17  and  18  show  that  the 
increase  in  frequency  due  to  pressure  is  approximately  the  same 
as  that  of  the  isotropic  case  for  9=0  and  90  degrees  (except 
for  mR/£  >  5  where  the  bending  stiffness  ratio  D{]/D22  begins  to 
take  effect). 

For  n  ;>  2,  the  frequencies  of  the  stiffened  cylinders 
under  internal  pressure  are  the  same  as  those  of  the  constant 
thickness  orthotropic  shell  (i.e.,  cylinders  1  and  5),  except 
for  large  values  of  xriR/l.  The  effect  of  pressure  on  the  two 
higher  frequencies  for  n  ^  1  is  negligible. 
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5.  CONCLUSIONS 


The  equations  for  the  vibration  of  a  thin  orthotropic 
cylindrical  shell  under  internal  pressure  were  deduced  from  the 
nonlinear  shell  theory  of  Washizu.  Application  was  made  to  four 
constant  thickness  orthotropic  cylinders,  two  stiffened  cylinders, 
and  for  purposes  of  comparison,  an  isotropic  cylinder.  This  in¬ 
vestigation  has  been  principally  concerned  with  the  effects  of 
circumferential  and  axial  stiffening,  as  well  as  internal  pres¬ 
sure,  on  the  frequencies  and  mode  shapes  of  an  orthotropic  cyl¬ 
inder  with  freely-supported  ends. 

The  significant  conclusions  are  that  there  is  a  difference 
in  behavior  of  the  frequencies  for  the  n  =  0,  n  =  1,  and  n  ^  2 
modes.  For  the  n  ^  2  modes,  the  lowest  frequency  is  associated 
with  predominantly  radial  motion.  For  the  n  =  0  mode,  the  lowest 
frequency  is  associated  with  either  the  uncoupled  torsional  mo¬ 
tion  or  the  predominately  radial  motion,  depending  on  the  axial 
wave  length.  For  the  n  *  1  mode,  there  is  strong  coupling  be¬ 
tween  circumferential  and  radial  motion,  and  the  lowest  fre¬ 
quency  is  associated  with  a  combination  of  circumferential  and 
radial  motion. 

For  the  n  ;>  2  modes,  the  frequency  of  the  predominantly 
radial  frequency  was  significantly  increased  by  circumferential 
stiffening,  while  this  frequency  was  unaffected  by  axial  stiffen¬ 
ing  for  all  but  large  values  of  mR/i.  For  the  n  »  0  mode,  there 
was  a  slight  increase  in  the  frequencies  of  the  torsional  mode 
and  predominantly  radial  modes  due  to  circumferential  stiffening, 
while  the  frequency  of  the  predominantly  axial  mode  was  unchanged. 
For  axial  stiffening,  the  frequency  of  the  predominantly  axial 
mode  was  significantly  increased,  while  the  frequency  cf  the  tor¬ 
sional  and  predominantly  radial  modes  was  slightly  affected.  For 
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the  n  *  1  mode,  circumferential  stiffening  increased  the  fre¬ 
quency  only  for  large  values  of  mR/l,  while  axial  stiffening  in¬ 
creased  the  frequency  mainly  for  small  mR/i . 

The  effect  of  pressure  on  the  frequencies  was  also  in¬ 
vestigated.  For  the  n  ^  2  modes,  for  both  circumferential  and 
axial  stiffening,  the  predominantly  radial  frequencies  were  sig¬ 
nificantly  increased  by  pressure.  For  the  n  ■  0  mode,  for  both 
circumferential  and  axial  stiffening,  the  torsional  frequency 
was  increased  by  pressure  at  small  mR/i  and  the  predominantly 
radial  frequency  was  increased  by  pressure  at  large  mR/i.  The 
predominantly  axial  frequency  was  not  affected  by  pressure.  For 
the  n  *  1  mode,  for  circumferential  stiffening,  the  frequency 
was  increased  by  pressure  at  small  mR/i,  while  for  axial  stiffen 
ing  the  frequency  was  slightly  increased  by  pressure. 

For  purposes  of  comparison,  the  frequencies  of  two 
stiffened  cylinders  were  also  investigated.  It  was  found  that 
the  frequencies  of  the  stiffened  cylinders  whose  bending  stiff¬ 
ness  ratio  ^'n/^22  was  ta^en  e<Iual  t0  stiffness  ratio 
bll//b22  of  the  constant  thickness  orthotropic  cylinder  are  essen 
tially  the  same  as  the  frequencies  of  the  isotropic  cylinder  for 
n  ■  0  and  n  ■  1  modes,  while  for  the  n  ^  2  modes  the  stiffened 
cylinder  frequencies  are  the  same  as  the  frequencies  of  the  con¬ 
stant  thickness  orthotropic  cylinder.  This  occurs  since  the 
n  *  0  and  1  modes  correspond  to  primarily  stretching  motion, 
while  the  n  ;>  2  modes  correspond  to  predominant ?y  bending  motion 
The  above  behavior  is  true  for  both  the  pressurized  as  well  as 
the  unpressurized  cases . 
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TABLE  1 

(1/R2) (D^ 2/^22)  Terms  for  Various  Thin  Cylindrical  Shell  Theories 
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TABLE  1  (Continued) 
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TABLE  2 

Initial  Stressing  Functions  for  Various  Thin  Cylindrical  Shell  Theories 
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TABLE  3 

Properties  of  Example  Constant  Thickness  Orthotropic  Cylinders 
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FIG.  2  EFFECT  OF  h/R  ON  FREQUENCY  DISTRIBUTION 
OF  n=0  MODE  FOR  ISOTROPIC  CYLINDER 
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FIG.  4  FREQUENCY  DISTRIBUTION  OF  n  =  0  MODE  FOR 
CIRCUMFERENTIAL  STIFFENING  -  CONSTANT 
THICKNESS  ORTHOTROPIC  CYLINDER 


48 


=  Rcu 


rrm/i 

FIG.  5  FREQUENCY  DISTRIBUTION  OF  n  =  0  MODE  FOR 
CIRCUMFERENTIAL  STIFFENING  -  CONSTANT 
THICKNESS  ORTHOTROPIC  CYLINDER 
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FIG.  7  FREQUENCY  DISTRIBUTION  OF  n  =  0  MODE  FOR 
AXIAL  STIFFENING -CONSTANT  THICKNESS 
ORTHOTROPIC  CYLINDER 


FREQUENCY  DISTRIBUTION  OF  n  =  0  MODE  FOR 
AXIAL  STIFFENING- CONSTANT  THICKNESS 
ORTHOTROPIC  CYLINDER 
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FIG.  9  FREQUENCIES  OF  n>|  MODES  FOR  CIRCUMFERENTIAL 
STIFFENING -CONSTANT  THICKNESS  ORTHOTROPIC  CYLINDER 
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FIG.  10  FREQUENCIES  OF  n>  I  MOOES  FOR  CIRCUMFERENTIAL 
STIFFENING -CONSTANT  THICKNESS  ORTHOTROPIC  CYLINDER 
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FIG.  12  FREQUENCIES  OF  n>|  MODES  FOR  AXIAL 
STIFFENING- CONSTANT  THICKNESS 
ORTHOTROPIC  CYLINDER 
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FIG.  13  FREQUENCES  OF  n>l  MODES  FOR  AXIAL 
STIFFENING -CONSTANT  THICKNESS 
ORTHOTROPIC  CYLINDER 
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FIG.  14  FREQUENCIES  OF  THE  BEAM  TYPE  MODE  (0=1) 
FOR  CIRCUMFERENTIAL  AND  AXIAL  STIFFENING 
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FIG.  15  FREQUENCY  DISTRIBUTION  OF  n  =  0 
MODE  FOR  CIRCUMFERENTIAL 
STIFFENING- STIFFENED  CYLINDER 
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FIG.  16  FREQUENCY  DISTRIBUTION  OF  n  =  0  MODE  FOR 
AXIAL  STIFFENING -STIFFENED  CYLINDER 
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FIG.  20  VARIATION  OF  FREQUENCY  WITH  PRESSURE -STIFFNESS 
PARAMETER  FOR  VARIOUS  CIRCUMFERENTIAL  WAVE 
NUMBERS  n-  ISOTROPIC  CYLINDER 
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FIG.  21  VARIATION  OF  FREQUENCY  WITH  PRESSURE  -  STIFFNESS 
PARAMETER  FOR  VARIOUS  CIRCUMFERENTIAL  WAVE 
NUMBERS  n- ISOTROPIC  CYLINDER 
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FIG.  22  EFFECT  OF  INTERNAL  PRESSURE  ON  FREQUENCIES 
OF  n  >  |  MOOES -CONSTANT  THICKNESS 
ORTHOTROPIC  CYLINDER 
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FIG.  23  EFFECT  OF  INTERNAL  PRESSURE  ON  FREQUENCIES 
OF  n>  |  MODES- CONSTANT  THICKNESS 
ORTHOTROPIC  CYLINDER 
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FIG.  24  EFFECT  OF  INTERNAL  PRESSURE  ON 
FREQUENCIES  OF  n  >  |  MODES  FOR 
ISOTROPIC  CYLINDER 
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FIG.  25  EFFECT  OF  INTERNAL  PRESSURE  ON  FREQUENCIES 
OF  n>  I  MOOES -CONSTANT  THICKNESS 
ORTHOTROPIC  CYLINDER 
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